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ABSTRACT 

The super-Earth GJ1214b transits a nearby M dwarf that exhibits 1% intrinsic variability in the 
near-infrared. Here, we analyze new observations to refine the physical properties of both the star and 
planet. We present three years of out-of-transit photometric monitoring of the stellar host GJ1214 from 
the MEarth Observatory and find the rotation period to be long, mostly likely an integer multiple of 53 
days, suggesting low levels of magnetic activity and an old age for the system. We show such variability 
will not pose significant problems to ongoing studies of the planet's atmosphere with transmission 
spectroscopy. We analyze 2 high-precision transit light curves from ESO's Very Large Telescope along 
with 7 others from the MEarth and FLWO 1.2 meter telescopes, finding physical parameters for the 
planet that are consistent with previous work. The VLT light curves show tentative evidence for spot 
occultations during transit. Using two years of MEarth light curves, we place limits on additional 
transiting planets around GJ1214 with periods out to the habitable zone of the system. We also 
improve upon the previous photographic U-band estimate for the star, finding V — 14.71 ± 0.03. 
Subject headings: eclipses — planetary systems: individual (GJ 1214b) — stars: individual (GJ 1214) 
— stars: low-mass, brown dwarfs — stars: rotation — stars: spots 



1. INTRODUCTION 

The transiting exoplanet GJ1214b offers an unparal- 
leled opportunity to explore the physical properties of 
super-Earth planets. With a mass (M p = 6.6 M®) and 
radius (R p = 2.7 R®) between those of Earth and Nep- 
tune, and a likely equilibrium temperature (T eq = 500K) 
cooler than for most transiting planets, GJ1214b repre- 
sents an intr iguing new kind of world w ith no Solar Sys- 
tem analog (;Char bonneau et al.l 120091) . Given intrinsic 
degeneracies in the mass-radius diagram in this regime 
dSeaeer et al.ll2007t lAdams et al.l[200l iRogers fc Seagerl 
l2010a| ). the bulk composition of the planet cannot be 
uniquely determined from current me asurements of the 
mass a nd radius alone. For example, IRogers fc Seagerl 
(20101)) can explain the observed mass and radius to 
within lcr with any of three generic physical models: 
(i) a mini-Neptune that accreted and maintained a low- 
mass H/He layer from the primordial nebula, (ii) a su- 
perfluid water-world with a sublimating H 2 envelope, 
or (iii) a rocky planet with an H-dominated atmosphere 
formed by recent outgassing. De tailed calculations 
of GJ 1214b's thermal evolution by Nettelman n et all 
(2010) favor a metal-enriched H/He/H20 envelope, find- 
ing that a water-only atmosphere would require an im- 
plausibly large water-to-rock ratio in the planet's inte- 
rior. 
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Fortunately, because GJ1214b transits a very nearby 
(13 pc), bright (K = 8.8), low-mass M dwarf (0.16 Mq), 
it is amenable to follow-up observations that could dis- 
tinguish among these hypotheses. In particular, the large 
(D = 1.4%) transit depth enables tran s it stu dies of the 
planet's atmosphere. I Miller- Ricci et all (|2009l ) show that 
measuring the amplitude of the planet's transmission 
spectrum (i.e., the wavelength-dependence of the tran- 
sit depth AD(X) caused by absorption at the limb of 
the planet) constrains the mean molecular weight of its 
atmosphere and, in turn, the hydrog en content of its 
outer e nvelope. Cases (i) or (iii) of IRogers fc Seagerl 
(|2010bfl would produce AD(X) w 0.1% variations in 
the transit depth across wavelengths accessible from 
the ground as well as Hubble and Spitzer Space Tele- 
scopes, while case (ii) 's dense atmosphere would result 
in variations below the sensiti vity of current instruments 
(jMiller-Ricci fc Fortnevll2l)10l ). 

Providing a potential complication, however, the host 
star GJ1214 shows roughly sinusoidal photometric mod- 
ulations that are presumably due to an asymmetric dis- 
tribution of spots on a rotating star. Such spots can 
bias planetary parameters as measured from transit light 
curv es whether or not they are occulted by the planet 
(e.g. lPont et~aLll2007t iDesert et al.ll2011bl ). partially de- 
coupling the observed transit depth D from the actual 
planct-to-star radius ratio R p /R+. Of particular impor- 
tance for transmission spectroscopy studies, the change 
in transit depth induced by spots can vary with both 
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time and wavelength, potentially mimicking the signal of 
a planetary atmosphere. 

Stellar spots have been observed in sever al transiting 
exopl anet systems around active stars (see iStrassmeierl 
l2009f) . Hubble S pace Telescope (HST) photometry 
(iRabus et al.1 12009ft a nd later ground-based follow-up 
(|Dittmann et all l2009h of TrES-lb has shown evidence 
for spot occultations in transit light curves. The high 
photometric precision and continuous coverage provided 
by the CoRoT satellite enabled detailed modeling of 
spotty transit and out- of-transit light curves for the hot 
Jupit ers CoRoT-2b fWolter et al.1 120091: iCzesla et al ' 
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2010) and CoRoT-4b (jAigrain et al 
2009). For the former, joint fits to 



the transit and out-of-transit flux showed that initial 
estimates of the pl anet's Rp/R* were 3% (9a) too low 
(jCzesla et al.1 12009) . The interpretation of the transit- 
ing super-Earth CoRoT-7b is obfuscated by the fact 
that both the transit depth and the reflex motion are 
well b elow the amplitude of activity- i nduce d modula- 
tions (jLeger et all 120091: iQueloz et ail I2009D . Reanal- 
yses of the CoRoT-7 radial veloc ities find changing 
values for the mass of CoRoT-7b (iHatzes et al.l 120101 : 
iFerraz-Mello et al.ll20Tol : lLanza et al.ll2010f ) and "call into 
question the significance of the mass measurements for 
both CoRoT-7ba nd the claimed outer planet CoRoT-7c 
(|Pont et al.ll20Tfj|) . 

Like GJ1214b, the well-studied hot Jupiter HD189733b 
(jBouchv et al.l 120051 ) is an ideal system for character- 
ization studies, but requires corrections for stellar ac- 
tivity. The host HD189733 is an active K2 dwarf 
(jMoutou et al.l I2007D with 2% peak-to-peak variabilit y 
in the optical (ICroll et al.ll2007t iMnler-Ricci et aT1l2008D . 
iHenry fc Winnl (|2008f ) undertook a long-term photomet- 
ric monitoring campaign from which they measured 
the 12 day stellar rotation period of HD189733. Ex- 
trapolation from their out-of-eclipse photometric spot 
characterization was useful for interpreting transmission 
spectroscopy results of individual t ransits from Hubble 
(jPont et al.ll2007t iSwain et al.l 120081 ) and measurements 
of the ther mal phase curve from Spitzer (K nutson et al.1 
120071 l2009T l . Understanding the time- variable surface of 
the star was even more crucial for broadband transmis- 
sion spectroscopy studies that rel y on comparing tra n- 
sit depths at dif ferent epochs (e.g.. iDesert et al.ll2009l : ?; 
iSing e t al. 2009); interpretation o f these data rely heavil y 
on the photometric monitoring of IHenry fc Winnl (|2008Q . 

To aid ongoing and future studies of GJ1214b, we 
present new data (J3) to characterize the star GJ1214's 
variability and estimate its rotation period (331). We com- 
pare the measured variability to a simultaneous analysis 
of 2 high-precision transit light curves from ESO's Very 
Large Telescope with 7 other new or previously published 
transits (j2J. Additionally, we place upper limits on the 
radii of other possible transiting planets in the system 
(SJS|) and present a rehned estimate of the star's V flux, 
which bears directly upon its metallicity as estimated 
using Mk and V — K relations. Finally, we discuss the 
implications of the measured variability for the proper- 
ties of the star and for transmission spectroscopy studies 
of GJ1214b's atmosphere (©. 

We also note the following correction. In 
iCharbonneau et"al"1 (|2009f) . we quoted a systemic radial 
velocity for GJ1214 that had a typo in the sign; the ac- 



tual velocity is 7 = +21.1 ± 1.0 km s 1 (i.e. a redshift). 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. MEarth Photometry 

We monitored the brightness of the GJ1214 system 
at a variety of cadences with the MEarth Observa- 
tory at Mt. Hopkins, AZ throughout the 2008, 2009, 
and 2010 spring observing seas ons. As described in 
iNutzman & Charbonneaul (|2008l ). the MEarth Observa- 
tory was designed to detect transiting exoplanets around 
nearby M dwarfs, and consists of eight identical 40- 
cm telescopes on German Equatorial mounts in a sin- 
gle enclosure at the Fred Lawrence Whipple Observa- 
tory (FLWO). Each telescope is equipped with a thinned, 
back-illuminated 2048x2048 CCD with a pixel scale of 
0.757" /pixel for a 26' held of view. For the bulk of the 
data presented in this work telescopes were equipped 
with a fixed, custom, 715 nm long-pass filter; the re- 
sponse is similar to a combination of the Sloan i + z band- 
passes and will be hereafter referred to as the "MEarth" 
bandpass. The MEarth Observatory is almost fully auto- 
mated and operates on every clear night, observing target 
stars selected from a list of 2000 nearby late M dwarfs 
(|Nutzman fc Charbonneaul [2008). In typical operating 
mode, each telescope observes its own list of 20-30 stars 
per night when they are above airmass 2 with the ca- 
dence and exposure times necessary to detect a transiting 
planet as small as 2 R^ in each star's habitable zone. 

Light curves are extracted automatically from MEarth 
images by a mod ified version of the Monitor pipeline 
(|Irwin et al.l 120071 ). using nightly flat field (dawn and 
dusk), dark, and bias exposures for calibration. A dif- 
ferential photometry correction for each frame was cal- 
culated from a robust, weighted fit to 78 automatically 
selected field comparison stars within 2.3 instrumental 
magnitudes of GJ1214. The mean MEarth - K color of 
these stars is 0.98; none were as red as GJ1214 (MEarth 
- K — 2.19). Predicted uncertainties for each measure- 
ment are calculated from a standard CCD noise model. 

GJ1214 was observed with three main cadences. 
"Low" cadence (20-30 minutes between exposures) was 
that associated with the normal survey mode, and was 
employed in the 2008 and 2009 seasons. "Medium" ca- 
dence (5-10 minutes) was implemented after the discov- 
ery of GJ1214b and was intended to boost sensitivity 
both to other transiting planets and to characterizing 
the out-of-transit variability of the star. "High" cadence 
(40 seconds) was employed at predicted times of transit 
to determine the system parameters. High cadence tran- 
sits were observed simultaneously with 7 or 8 MEarth 
telescopes for greater precision, as the systematic noise 
sources and scintillation patterns among pairs of MEarth 
telescopes appear to be largely uncorrelated, so the S/N 
improvement scales with the square root of the number 
of telescopes. In Table [TJ we present one new MEarth 
transit light curve, along with the four MEarth and two 
KeplerCam light curves t hat were analyzed bu t not m ade 
electronically available in ICharbonneau et all (|2009| ) . 

While we include data from 2008 for the rotation anal- 
ysis, we caution that these observations took place dur- 
ing MEarth's early commissioning, before the observ- 
ing strategy and software were finalized. Changes to 
the telescope throughout the season may have corrupted 
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the season-long stability. Importantly, a field acquisition 
loop designed to mitigate fiat-fielding errors by bringing 
each star back to the same pixel was not implemented 
until the late spring of 2008. During the 2008 season, 
a Bessell-prescription I filter (Besseli Hl990| ) was used in- 
stead of the custom MEarth bandpass. 

2.2. V-band Photometry from KeplerCam. 

An independent out-of-transit light curve was obtained 
through Harris V and I filter^ with KeplerCam on the 
1.2 meter reflector at FLWO atop Mt. Hopkins, AZ. 
The observations discussed here were gathered in service 
mode from 26 March 2010 until 17 June 2010, after 
which date the mirror was taken off and put back on the 
telescope, introducing an uncorrectable systematic offset 
to the field light curves so later data had to be discarded. 

Given the large night-to-night positional shifts of the 
field, we made an effort to quantify and ameliorate flat- 
fielding errors by sampling multiple regions of the de- 
tector, with each observation consisting of a set of three 
exposures offset by 3 arcminute dithers. Individual ex- 
posures had theoretical noise limits ranging from 0.3% 
to 2%, but the scatter among dither points suggested 
that calibration errors from flat-fielding introduced a 1% 
noise floor to the light curve. Dark sky flats were gener- 
ated and corrected over time for changes at high spatial 
frequencies (i.e. dust donuts) by nightly dome flat expo- 
sures. 

We measured calibrated V and / magnitudes (Table [5]) 
to imp rove on previously pu blished photograp hic esti- 
mate s (jLepine k. Sharal [20051 ). Standard fields ()Landoltl 
1992) were observed on the nights of 26, 27, and 28 
March. Conditions were clear, although seeing as poor as 
10" FWHM was witnessed. We estimate the calibration 
uncertainties for the nights from the scatter in multiple 
standard exposures. 

2.3. Light Curves from VLT-FORS2 

Spectra of GJ1214 and 6 comparison stars were gath- 
ered during three transits of GJ1214b using ESO Direc- 
tor's Discretionary Time on the VLT (Prog. ID #284.C- 
5042 and 285.C-5019). As described by iBean et al.l 
(pOioi) . the primary purpose for obtaining these data 
was to measure the transmission spectrum of GJ1214b's 
atmosphere by generating multiwavelength transit light 
curves and determining the wavelength-dependence of 
the transit depth. In this work, we generate and analyze 
high precision "white" light curves by summing together 
all the photons collected in each spectrum. 

Observations were performed in queue mode with 
the multiob.ject, low d ispersion spectrograph FORS2 
(jAppenzeller et al.l 119981) on VLT/UT1. The spectro- 
graph was configured with the 600z+23 grism with a cen- 
tral wavelength of 900 nm and the red-sensitive (MIT) 
CCD in the standard 2x2 read mode. Exposure times 
were 20-40 seconds, and the readout time was 37 seconds. 
A custom slit mask was used; each slit was a rectangle 
12" in the dispersion direction and 15-30" in the spatial 
direction, small enough to isolate GJ1214 and the com- 
parison stars but large enough that changing slit losses 
due to variable seeing were negligible. Wavelength cali- 
bration exposures with a He, Ne, Ar emission lamp were 

2 http://www.sao.arizona.edu/FLWO/48/CCD.filters.html. 



taken through a 1" slit the day after each set of observa- 
tions. Given the position of the comparison stars on the 
chip, the wavelength range 780 to 1000 nm was used in 
this analysis. The CCD response governs the red edge of 
this range, and the spectral response is similar to that of 
the MEarth bandpass. 

After bias subtraction and flat fielding, we used the 
comparison stars to correct for the time varying zero- 
point of the system. For each exposure, we extracted 
ID sp ectra from th e images using the optimal extrac- 
tion of Home (1986), and divided the total flux (summed 
over wavelengths) of GJ1214 by the total flux in all the 
comparison stars. Theoretical error bars calculated from 
photon statistics alone were assigned to each point. Each 
exposure yielded 1 — 3 x 10 8 photons from GJ1214 and 
twice as many from comparison stars. 

The corrected GJ1214 light curves exhibit systematic 
trends which we correct for by fitting a second-order 
polynomial function of time. In 2] we propagate the 
uncertainty from the systematics corrections through to 
the transit parameters. We searched for correlations be- 
tween the relative flux and airmass, seeing, and positional 
shifts in the dispersion and cross-dispersion directions. 
The relationships were more complicated than low-order 
polynomials, so we did not attempt to rem ove them us- 
ing c ommon decorrelation techniques (e.g. I Burke et al.l 
1201(3) . 

We also tested whether the observed drifts in flux could 
be caused by a changing color-dependence of the atmo- 
spheric extinction along the line of sight. To do so, we 
applied differential photometry corrections to individual 
spectral channels before combining them, to allow each 
wavelength its own extinction. This procedure did not 
remove the systematic trends. 

We suggest the follow ing as a more probab le explana- 
tion for the systematics. IMoehler et al.l (|2010f ) found that 
the linear atmospheric dispersion corrector (LADC) on 
the telescope has surface features that affect is sensitivity 
across the field of view. Because the LADC is positioned 
before the field rotator in the optical path and rotates 
relative to the sky, individual stars can drift across these 
features and encounter throughput variations that are 
not seen by the other comparison stars. No rotationally- 
dependent flat - fields were applied to these data, although 
IMoehler et al.1 (|2010[) provide a route to a possible cor- 
rection. 

The first two "white" light curves, normalized to their 
median out-of-transit flux level, are published in Table Q] 
and shown in Fig. SJ 

A third transit observation was attempted on 2010 Jul 
22. The brightest comparison star could not be used be- 
cause it saturated mid-transit. The exposure times were 
cut in half immediately after egress, and a notable off- 
set is visible in the transit light curve, perhaps due to 
an uncorrected non-linearity in the detector. Although 
the light curves of the first two transits were robust to 
the choice of comparison stars, the third changed signif- 
icantly depending on which set of comparison stars was 
used. We show this transit in Fig. El but exclude it from 
all following analyses. 

3. ROTATION PERIOD OF GJ1214 

With a growing understanding of the systematic effects 
present in MEarth data, we revisit the issue of GJ1214's 
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TABLE l 
Transit Light Curves 



Mid-exposure Time a Relative Flux b Error c Airmass Instrument Cycle d 



2454964.8926699 


0.99622 


0.00317 


1.1197 


MEARTH 





2454964.8933879 


0.99777 


0.00321 


1.1198 


MEARTH 





2454964.8941049 


1.00578 


0.00325 


1.1199 


MEARTH 





2454980.7148766 


0.99733 


0.00185 


1.6190 


FLWO 


9 


2454980.7153966 


1.00044 


0.00185 


1.6140 


FLWO 


9 


2454980.7158946 


0.99959 


0.00185 


1.6080 


FLWO 


9 


2455315.7660750 


0.99786 


0.00014 


1.2125 


VLT 


221 


2455315.7668492 


0.99762 


0.00014 


1.2105 


VLT 


221 


2455315.7702432 


0.99763 


0.00012 


1.2020 


VLT 


221 



Note. — This table is presented in its entirety in the electronic edi- 
tion; a portion is shown here for guidance regarding its form and content. 
a Times are given as BJDtd r, Barycentric Julian Dates in the Barycen- 
tric Dynamical Time system (jEastman et al.|[2010f) . 
b Differential photometry corrections have been applied, but additional 
systematic corrections (see text) have not. Each light curve has been 
divided by the median out-of-transit flux. 

c Theoretical la errors have been calculated from a standard CCD noise 
model. 

d Time measured from the reference epoch in units of the orbital period. 



TABLE 2 
Photometry of G J 1214 



Filter 


Magnitude 


Source 


V 


14.71 ± 0.03 


this work 


1 


11.52 ±0.03 


this work 


J 


9.750 ± 0.024 


2MASS a 


H 


9.094 ±0.024 


2MASS a 


K 


8.782 ± 0.024 


2MASS a 



a ISkrutskie et all (|2006D 

intrinsic variability. In th e discovery paper for GJI2I4b 
(jCharbonneau et al.|[2009h . wc stated that the dominant 
periodicity seen in the out-of-transit light curve of the 
star GJ1214 had an 83 day period, implying that this 
was the rotation period of the star. Here, we revisit the 
question of GJ1214's rotation period with another season 
of observations. 

Semi-stable spot complexes on the surface of a star im- 
print photometric modulations that can be approximated 
as a sinusoid with a fundamental period that matches 
the stellar rotation period. We search each year's light 
curve with a weighted, least-squares periodogram that 
has been modified to simultaneously fit for stellar vari- 
ability along with scaled templates of known systematic 
effects. These systematics are discussed in the next sev- 
eral subsections. To account for the likely evolution of 
spots with time, we investigate the 2008, 2009, 2010 data 
sets separately and do not require a coherent sine curve 
to persist over multiple years' data. 

3.1. Avoiding Persistence in MEarth Light Curves 

The MEarth detectors are subject to image persis- 
tence; pixels that are illuminated in one exposure can 
show enhanced dark current in subsequent exposures, 
which decays exponentially with a half hour time scale. 
Because the dark current in a given pixel depends on how 
recently that pixel was illuminated, differential photom- 



etry light curves can show baseline shifts between obser- 
vations taken at different cadences, as well as 'ramps' at 
the start of a high-cadence sequence of exposures. 

Correcting for these changing baseline shifts would re- 
quire a simultaneous modeling of the complete photon 
detection history of every pixel and is impractical. Dur- 
ing MEarth's normal survey mode, we purposely center 
subsequent targets on different pixels to avoid persistent 
charge stacking up. As the effect is most noticeable for 
data with cadence shorter than 5 minutes, we circumvent 
the problem by throwing out from the rotation period 
analysis all but the first point of any segment with such 
cadence. 

3.2. Adding a Systematic 'Jitter 1 

Even having removed the highest cadence data, the 
sampling in the MEarth GJ1214 light curve can vary 
from N=l to N=35 points per night, with a typical the- 
oretical noise limit per point of 3 millimagnitudes. In a 
strict least-squares sense, if the noise in our data were 
accurately described by an uncorrelated Gaussian pro- 
cess (e.g. photon noise), then every exposure should be 
allowed to contribute on its own to the period search, 
meaning the uncertainty associated with each night 
would go down as 1/y/W. Although MEarth telescopes 
have achieved such white noise d own to millimagni- 
tude levels within individual nights dCharbonneau et al.l 
I2009T ) . photometric variations between nights are most 
likely dominated by subtle changes in the telescope that 
are not corrected by our calibration efforts. The 1/y/N 
weighting scheme would unfairly bias a period search to 
fit only a few well-sampled nights. 

To account for this in each light curve, we remove all 
in-transit exposures and bin the data to a nightly time 
scale. For each night with N data points, we calculate 
an inverse variance weighted mean flux and time, using 
the theoretical errors calculated for each exposure in the 
weighting. To each of these nightly bins we assign an 
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Fig. 1. — Periodograms showing the A% 2 improvement achieved 
by fitting a sinusoid of a given period + common mode + meridian 
flip over a null model consisting only of common mode + meridian 
flip for MEarth (top 3 panels) and FLWO V-band photometric 
monitoring (bottom). Possible periods as short as 1 day (left) and 
a zoom in to longer rotational periods (right) are shown. Periods 
for which less than one full cycle is observed per season are shaded. 
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Fig. 2.— Out-of-transit light curves of GJ1214 from the MEarth 
Observatory (top three panels) and in V from the FLWO 1.2m (bot- 
tom panel). Individual exposures (gray points) and nightly binned 
values with errors that include the systematic 'jitter' (black cir- 
cles, see text for details) are shown for each. A sine curve at the 
proposed 53 day period (derived from 2010 MEarth data) is shown 
(solid lines) at the best-fit phase and semiamplitude for each light 
curve. The MEarth points have been corrected for the best-fit 
common mode and meridian flip decorrclation. 



a 



nightly 



bin 



a 



jitter 
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where Ohm is the intrinsi c stand ard error on the mean of 
the nightly bin (RMS/y/N — 1) and a jitter is a constant 
noise floor term to capture the night-to-night calibration 
uncertainty. 

The signal lost in the binning process should be min- 
imal. Preliminary searches of unbinned data and visual 
inspection of high cadence nights revealed no significant 
periodic signal at periods shorter than f day. Under 
the (untested) assumption that the stellar spin is roughly 
aligned with the orbital angular momentum, the upper 
limit on the projected rotation velocity of vsini < 2 km 
s _1 would correspond to a rotation period P ro t > 5 days. 

We use observations on successive nights to estimate 
a juter for the MEarth and the V-band light curves for 
each observing season. This assumes, on the basis of 
the apparent lack of short term variability, that the flux 
difference between pairs of nights is dominated by sys- 
tematics. We find a juter = 0.0052, 0.0058, 0.0038, 0.0067 
magnitudes for the 2008, 2009, 2010 MEarth and 2010 V 
light curves. As these values are comparable to the pre- 
dicted noise for most exposures, the quadrature addition 
of a juter means we weight most nights roughly equally. 

3.3. Correcting for Meridian Flips in MEarth Light 

Curves 

MEarth's German Equatorial mounts rotate the detec- 
tors 180° relative to the sky when switching from nega- 
tive to positive hour angles. Thus, the target and com- 
parison stars sample two different regions of the detec- 
tor. Given imperfect flat field corrections, an offset is ap- 
parent in many MEarth light curves between exposures 
taken on either side of the meridian. 

To account for this effect, we allow different sides of 
the meridian to have different zcropoints. We construct 
a meridian flip template m(ti), which for an unbinned 
light curve would consist of binary values corresponding 
to the side of the meridian at each time stamp By 
extension, for each night of the binned data U, we define 



m(ti 



7(n H 



(2) 



where n + and rt_ are the number of data points with 
positive or negative hour angles in a given nightly bin. 
We allow a scaled version of this template to be fit si- 
multaneously with the period search. 

3.4. Correcting for Water Vapor in MEarth Light 
Curves 

Because the wide MEarth bandpass overlaps signifi- 
cant water absorption features in the telluric spectrum, 
the color-dependence of the throughput of our observ- 
ing system is sensitive to the precipitable water vapor 
(PWV) in the column overhead. The fraction of stellar 
photons lost to water vapor absorption from a typical 
MEarth target M dwarf is much larger than the frac- 
tion lost from the (typically solar-type) comparison stars. 
When PWV along the line of sight to a star varies, a cru- 
cial assumption of simple differential photometry - that 
stars are experiencing the same losses - is violated. 

Although the variations in any particular light curve 
might come either from the PWV induced noise or from 
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intrinsic stellar variability, we can harness the ensem- 
ble of M dwarfs observed by MEarth at any particu- 
lar time to characterize and correct for this effect. To 
do so, we construct a "common mode" template by ro- 
bustly (median) binning all the differential photometry 
light curves of all M dwarfs observed on all eight MEarth 
telescopes into half hour bins. This averages out uncor- 
rected stellar variability and serves as an estimate of the 
atmospheric variation that is common to all red stars ob- 
served at a given time. We only use data during times 
when we have > 50 and > 30 targets contributing to 
a bin. The strongest periodicities in the common mode 
templates are 2 5.1 days for sprin g 2009 and 14.5 days 
for spring 2010 ([Irwin et al.ll2010l) . If left uncorrected, 
such periodicities could appear as spurious intrinsic stel- 
lar variability. 

To correct the GJ1214 light curve for this effect, we 
interpolate the "common mode" to the unbinned time 
stamps. We then perform the nightly binning on it to 
construct a common mode template c{U), which we use 
for simultaneous decorrelation (see next section). This 
binning is justified because the typical common mode 
variation within a night is typically at the level of 1-2 
millimagnitudes, much smaller than the night-to-night 
or week-to-week changes we hope to correct. 

3.5. Periodograms 

We generate a periodogram by calculating the \ 2 of a 
weighted linear fit of the light curve AF(ti) (in magni- 
tudes) to a model 

AF slne {U) = A sin ( 27T (ti-to) \ +B m ^ + c C ( U ) + D 

\ Prot / 

where A is the semiamplitude of the sinusoidal variabil- 
ity, to is an epoch, P rot is the stellar rotation period, B 
and C are scale factors for the systematics, and D is a 
constant offset. Wc compare this \ 2 to that of the null 
hypothesis that AF is explained by the systematics alone 

AF nuU (ti) — B m(ti) + C c(U) +D. (4) 

Mathematically, this procedure woul d be identical to tra- 
dition al least-squares periodogra ms (lLom b 1976; S cargld 
fl982f) if we fixed B = C = 0. In llrwin et all (|2010D . we 
use a similar method to estimate photometric rotation 
periods for a sample of 41 MEarth M dwarf targets and 
test its sensitivity with simulations. 

In Fig. [1] we plot the \ 2 improvement ( A% 2 = Xnu\\ ~ 
xline) between these two hypotheses for each of the three 
MEarth seasons and the short 2010 V-band campaign. 
Periods for which less than one cycle would be visible 
have been masked. 

The most prominent peak among all the periodograms 
in Fig. Q] is that at 53 days from the 2010 MEarth data, 
which corresponds to a semiamplitude of A = 3.5 ± 0.7 
millimagnitudes, where the uncertainty has been esti- 
mated from the covariance matrix of the linear fit. We 
estimate the false alarm probability (FAP) for this pe- 
riod by running the complete period search on 10 4 time 
series that consist of the best-fit scaled versions of m(ti) 
and c(U) and randomly generated Gaussian noise set by 
^nightly, recording the A\ 2 of the best peak from each 
iteration. We find that a FAP of 10~ 4 corresponds to 




0.010 0.005 0.000 -0.005 -0.010 
AMEarth 



Fig. 3. — For a sub-sample of nights measured in both MEarth 
and V bandpasses, the nightly bins plotted against each other with 
the assigned u\y in error bars. The best-fit slope (black line) and la 
interval (shaded region) are shown. 

A\ 2 — 28, much less than the achieved A\ 2 — 41. 

There is a nearby, but statistically insignificant, peak 
at 51 days in the 2009 MEarth data. Both 2008 and 2009 
MEarth light curves are dominated by long-period trends 
that are unresolvable in each year. Of the resolved peaks 
in the 2008, the strongest is at 81 days (FAP < 10" 4 ). 
In spite of the formal significance of this last peak, we 
caution that mid-season changes to the then still uncom- 
missioned observatory might also account for the vari- 
ations seen. One conclusion is robust; our 3 years of 
MEarth light curves show no evidence for any rotational 
modulation with a period shorter than 25 days. 

Given that 2010 had the most uniform sampling and 
cadence, we tentatively suggest P ro t = 53 days as our 
current best estimate of GJ1214's likely rotation period. 
Fig. [2] shows the 3 MEarth and 1 V binned light curves 
with a sinusoid whose period has been fixed to our esti- 
mated P rot = 53 days but whose amplitude and phase 
have been fitted to the data. The fit is acceptable for 
2009 MEarth, but clearly fails for the 2008 MEarth data. 

We stress the caveat that the true rotation period could 
instead be a longer multiple of our quoted 53 day period 
(e.g. P rot ~ 100 days) if the star exhibits multiple, well- 
spaced active regions. This kind of harmonic confusion 
appeared and was addressed in studies of Pr oxima Cen 
(jBenedict et al.lll998t iKiraga fc Stepienll2007f ). Prelimi- 
nary data for GJ1214 collected in 2011 while this paper 
was under review do not seem to show evidence for a 53 
day period, preferring instead a much longer one. Due to 
this factor-of-n uncertainty in the true rotation period, 
we do not quote a formal error bar on our 53 day period 
estimate. 

3.6. Chromatic Spot Variation 

By itself, the V-band light curve prefers a period of 
41 days (Fig. [2 FAP = 5xl0~4). When forced to fit a 
53 day period (Fig. [2]) , these data show a phase offset of 
only 2 days relative to the simultaneous MEarth data. 
The semiamplitude of this fit is 7 ± 3 millimagnitudes, 
twice that seen in the MEarth bandpass. 

There are 28 nights when observations were obtained 
in both MEarth and V band. In Fig. [3] we plot the nightly 
bins against each other; the apparent correlation suggests 
that the two instruments are observing the same stellar 
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TABLE 3 
Summary of Transit Light Curves 



Cycle 


UT Date 


Inst runic nt 


RMS (ppm) 


CeidcncG (sgc.J 





2009 May 13 


MEarth 


2950 


60 


9 


2009 May 29 


FLWO 


1960 


45 


9 


2009 May 29 


MEarth 


1580 


(binned) 45 


11 


2009 Jun 01 


FLWO 


2060 


45 


11 


2009 Jun 01 


MEarth 


1240 


(binned) 45 


21 


2009 Jun 17 


MEarth 


1620 


(binned) 45 


221 


2010 Apr 29 


VLT 


380 


72 


228 


2010 May 10 


MEarth 


1770 


(binned) 45 


233 


2010 May 18 


VLT 


350 


72 



variability across two bands and not telescope systemat- 
ics. We fit a line to the r elation, acco unting for errors in 
both AMEarth and AV (|Pressll2002D and find a slope of 
2.4±0.8. While the significance is marginal, we take this 
as further evidence that the amplitude of the variability 
in V is greater than that in the MEarth band. 

If starspots have a temperature (T.) that is only mod- 
estly lower than the stellar effective temperature (T c g), 
the color-dependence will arise from the spectrum of 
the spot rotating in and out of view. The factor of 2 
we see would be consistent with T,/T e g 90 — 95% 
as is commonly assumed in M dwarf ecl ipsing binaries 
(e.g. lMorales et alJl2009Ulrwin et al.ll2009f ). Totally dark 
spots (T. = OK) would produce less of a chromatic varia- 
tion, but would still be sensitive t o the spectral signat ure 
of the stellar limb-darkening fe.g. lPoe fc Eatonlll985n . 

4. FITTING THE TRANSIT LIGHT CURVES 

We perform a simultaneous fit to the 4 transit light 
cu rves from MEarth a n d 2 f rom KeplerCam published 
bv iCharbonneau et al.l fj2009f ) . 1 more transit collected 
by MEarth in spring 2010, and 2 high-precision transits 
from the VLT. This totals 9 light curves of 7 indepen- 
dent transits, as shown in Figures |3] and [5] We employ 
a model corresponding to a circular planet transit ing a 
smooth, limb-darkened star (|Mandel fc Agoll 120021 ) that 
has the following parameters: the planct-to-star radius 
ratio R p /R+, the stellar radius i?*, the total transit dura- 
tion <i4, two quadratic stellar limb darkening parameters 
Mi and U2 for each of the 3 telescope systems used, and 
7 mid-transit times T c . The reparameterization of the 
scaled semimajor axis a/i?* and inclination i in terms 
of R+ and £14 substantially reduces the degeneracies in 
the problem, lea ding to an efficient exploration of th e 
parameter space dBurke et al.[|2007l : iCarter et al.l l2008) . 

We fix the orbital period to P = 1.5804043 days 
(|Sada et al.l [2010h . Given the upper limit on eccen- 
tricity from radial velocities and the short circulariza- 
tion time for GJ1214b (10 6 years f or Q' p = 100 and 

Q'„ = 10 6 , following iRavmond et all 120081 ) . we assume 
an eccentricity e = throughout. Where necessary to 
derive physical parameters from the geometric parame- 
ters in the light curve fit, we adopt the th e stellar mass 
M* = 0.157 M (jCharbonneau et al.ll2009h : wc describe 
how we propagate the 0.019 M Q uncertainty on this 
value to the other errors in the next section. 

To account for the systematic trends present in the 
VLT light curves, we introduce a correction to the base- 
line stellar flux parameterized as a parabola in time 
(a + bt + ct 2 ) to each transit. Most of the MEarth 



and FLWO light curves showed no significant systematic 
trends, so only a single out-of-transit baseline flux level 
was fit to each night. The MEarth light curve on the 
night of 2010 May 10 showed a strong correlation with 
airmass, so we also included a linear trend with airmass 
for this one night. 

4.1. x 2 minimization 

We determine the best-fit values of the 30 model pa- 
rameters by using an implementation of the Levenber, 
Marquadt (LM) routine called MPFIT (jMarkward 
2009) to minimize the value of 



X 



N 



(-^S -^model) 



(5) 



where Fi are the N = 1495 flux measurements (Table Q}, 
o~i are their uncertainties, -Fmodoi is the model described 
above. 

After this initial fit, the uncertainty estimates for each 
light curve were increased until the reduced \ 2 01 the out- 
of-transit residuals was unity and the fit was repeated. 

While the LM fit provides a linearized estimate of 
the covariance matrix and errors of the parameters, 
this estimate is too precise because the method a) does 
not fully sample along non-linear correlations between 
highly-degenerate parameters and b) does not account 
for correlations between data points. In what follows, we 
calculate more conservative and realistic errors through 
a bootstrap method that addresses both these issues. 

4.2. Error Estimates by Residual Permutation 

Although the noise in the MEarth and FLWO light 
curves is well described by a white Gaussian process, 
the much lower photon noise in the VLT light curves re- 
veals underlying low level serial correlations among data 
points, or 'red' noise. The autocorrelation function of the 
VLT residuals is above 0.25 out to 4 and 2 data point lags 
respectively for the two transits, and binning the resid- 
uals by N points reduces the scatter more slowly than 
1/\^N. To quantify the excess uncertainty in the fitted 
parameters due to this red noise, w e perform a 'residual 
permutation ' boot s trap simulation dMoutou et al.l 120041 : 
iGillon et all l2007t iDesert et all |2009j) that fits resam- 
pled data while preserving the correlations between data 
points. This analysis is carried out simultaneously for 
the VLT, FLWO and MEarth light curves. 

After subtracting the best-fit model from the ensem- 
ble of light curves, we perform 2 x 10 4 iterations of the 
following procedure. Preserving the time stamps for all 
the exposures, we cyclicly permute the residuals for each 
light curve by a random integer (shifting along the series 
and wrapping back from the last exposure to the first), 
inject the best-fit model back into the set of shifted resid- 
uals, perform the LM fit on the simulated light curve 
ensemble, and record the results. To initialize the pa- 
rameters for the LM fit, we select guesses drawn from 
a multivariate Gaussian distribution whose covariance 
matrix is scaled up by 2 2 (i.e. 2a) from the LM's esti- 
mate. We excised those bootstrap samples that found 
a best fit with unphysical limb- darkening paramet ers 
(ui < 0, tii + u 2 > 1, ui + 2u 2 < 0, iBurke et ai1l2007l ). 

In Table |U we quote the best-fit value for each pa- 
rameter and uncertainty bars that exclude the lower and 




Fig. 4. — Transit light curves from summed FORS spectra from the VLT, before (top) and after (middle) subtracting a second-order 
polynomial function of time. The third transit (open circles) exhibited larger uncorrectable systematics than the first two (see text), and 
was excluded from further analysis. Residuals from the quoted model (gray lines) are also shown (bottom), with dotted lines indicating the 
in-transit duration. 



lot . 




-50 -25 25 50 -50 -25 25 50 -50 -25 25 50 -50 -25 25 50 -50 -25 25 50 -50 -25 25 50 -50 -25 25 50 

Time from Mid-Tronsit (minutes) 



Fig. 5. — Transit light curves from MEarth and KeplerCam. Residuals from the quoted model (gray lines) are also shown (bottom), with 
dotted lines indicating the in-transit duration. The cycle 228 MEarth transit shows a linear trend with airmass that is included in the fit. 



upper 15.9% of the bootstrap samples (i.e. the central 
68.3% confidence interval), where parameters that were 
not directly fit have been calculated analytical ly from 
those that were (|Seager fc Mallen-Ornelasll200l . Hav- 
ing fit quadratic limb-darkening parameters for each in- 
strument, we present similar confidence intervals of the 
coefficients u\ and U2 in Table [SJ 

Fig. [6] summarizes the correlations among the param- 
eters R p /R±, a/Ri,, i, and the linear combination of the 
limb-darkening parameters u\ + ui for the three tele- 
scope systems. As the difference between central and 
limb specific intensities, u\ + ui correlates strongly with 
the transit depth and R p /R+. We also show in Fig [5] the 



bootstrap histograms for these parameters. 

Although we use i?* as a fit parameter, the quantity 
that is actually constrained by the light curves is p*. 
To calculate the true uncertainty on R+ (and R p ), we 
calculate for all of our bootstrap samples using the 
fixed M+ = 0.157M Q , assign values of M* drawn from 
the appropriate Gaussian distribution, and recalculate 
R± from and M*. 

We find results consistent with ICharbonneau et "all 
(2009). Despite the high precision of the VLT light 
curves, the uncertainties for most parameters are compa- 
rable to the earlier work, and that on R p /R+ is slightly 
larger. This is due in part to the correlated noise anal- 
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# 
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0.114 0.116 0.118 14.2 14.6 15.0 15.4 88.5 89.0 89.5 

R„/R. a/R. i 



0.4 0.6 0.8 1.0 0.4 0.6 0.8 1.0 0.4 0.6 0.8 1.0 
NMEarth / . , \FLW0 / , , \VLT 



(u, + u 2 ) 



( Ul +u 2 r 



{ i ^v 
(u, + u 2 ) 



Fig. 6. — Matrix showing the correlations between selected transit parameters from a simultaneous residual permutation analysis of 5 
MEarth, 2 FLWO, and 2 VLT light curves. Histograms are shown for each parameter, as well as contours that contain 68% (dark gray) and 
95% (light gray) of the samples for each pair of parameters. For each telescope system, we show a linear combination of the limb-darkening 
coefficients ui + 112 that is strongly correlated with the other parameters. 



ysis we perform that Charbonncau ct al 



(200l) did not. 
More significantly, iCharbonneau et al.l (|_2009) fixed the 
quadratic limb-darkening parameters to theoretical val- 
ues while we fit for them directly. The extra degrees 
of freedom allowed by our relaxation of astrophysical 
assump tions are known to increase the un certainty on 
Rp/R* dBurke et al.ll2007tlSouthworthl l2008l). For com- 
parison, [ChaAonnea^^LaO~l200^ - used coefficients ap- 
propriate for a 3000K, log g = 5, PH OENIX atmo sphere, 
specifically Cousins / ("tabulate d in iBesselll Il990h coeffi- 
cients (iti = 0.303, u 2 = 0.561; IClaretl U998lT as an ap- 
proximation for the MEarth band pass and Slo an z co- 
efficients (u x = 0.114, it 2 = 0.693: IClaretl 120041 ) for the 
KeplerCam FLWO data. While our individual fitted val- 
ues differ from these, the integral over the stellar disk 
(I — ui/3 — U2/6) is very well reproduced. 

4.3. Transit Timing Results 



Mid-transit times (equivalent to times of inferior con- 
junction given the assumed circular orbit) are printed in 
Table [6] Our uncertainty estimate on the VLT transit 
times is 50% larger than it would be if we had ignored 
correlations in the data, but still only 6 seconds. This 
uncertainty is within a factor of two of the uncertainty on 
the highest-precision transit times yet measured from ei- 
ther the ground or space: the 3-second me asurements of 
HDI89733b with Spitzer (|Agol et al.ll2010l) . For stars of 
comparable brightness to GJf 214, large aperture ground 
based telescopes offer a powerful tool for precision transit 
times. 

A revised linear ephemeris is shown in Table U de- 
rived with a weighted leas t squares method from our 
transit times and those of iSada et al.l ([20 10). Residu- 
als are shown in Fig. [7] The reduced x 2 of the linear fit 
= 0.62 gives no indication of transit timing variations 
over the time scales probed, and a Bayesian model com- 
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TABLE 4 

GJ1214B Parameters and Uncertainties 



Parameter 



Previous Work 



This Work 3 



Rp j R* 
a/R* 
i (deg.) 
b 

P*(g/cm 3 ) 
R*(R©) b 

ti2 = t34 (minutes) c 

t23 (minutes) c 

tl4 (minutcs) c 

P (days) 

T (BJDtdb) 



0.1162 ± 0.0007 
14.66 ± 0.41 

88.62l°j| 

V.d04_ Q og2 

23.9 ±2.1 
0.2110 ± 0.0097 
2.68 ±0.13 



1.5804043 ± 0.0000005 
2454966.52506 ± 0.00006 d 



0.1171 ±0.0010 
14.93 ± 0.24 

88.80±°;»S 

n Qio+0.046 

25.2 ± 1.2 
2064+ - 0086 

U.ZUDt_ oog6 

2.64 ±0.13 

6.00ig;iS 
40.l0lg;« 

52 11+ - 25 
1.58040490 ± 0.00000033 
2454966.525042 ± 0.000065 



References. — ICharbonneau et all (|2OO90 . ISada et al.l (120101 ) 
a Best-fit value and confidence intervals that exclude the upper and 
lower 15.9% of the residual permutation bootstrap samples. Where 
upper and lower limits differ by > 10%, asymmetric errors are shown. 
b Uncertainty calculated as described in the text assuming M* = 
0.157±0.019M Q 

c Times between contact are tyi — £34 = ingress/egress time, t 2 3 
rl'irat ion of full eclipse, tu = total duration of transit, tu was fit 
directly from the light curve, and the others were calcuated from the 
fitted parameters (see iSeager fc Mallen-Ornelasl 120031 iCarter et al.l 
120011 . 

d Converted from BJD UTC to BJD TDB and added 1 cycle for ease 
of comparison. 



TABLE 5 

Inferred GJ1214 Limb-darkening Coefficients 



TABLE 6 
G J 12 14 Transit Times 



Instrument A range (nm) 


«i a 


«2 a 


MEarth 715-1000 
VLT-FORS2 780-1000 
KeplerCam z 850-1000 


0.53 ±0.13 
0.34 ±0.31 
0.26 ±0.11 


-0.08 ±0.21 
0.28 ±0.46 
0.26 ±0.19 



a Best-fit value and confidence intervals that ex- 
clude the upper and lower 15.9% of the residual 
permutation bootstrap samples. 



parison test (see lBurke et al.ll2010( ) does not show signifi- 
cant evidence for a model with linear and hypothetical si- 
nusoidal components over a purely linear ephemeris. We 
note that a hypothetical 1.0 M s planet in a 2:1 mean- 
motion resonance with GJ1214b would introduce ~ 100 
second transit timing variation s that woul d be easily de- 
tected in these data (following, Bean 2009|). 

4.4. Occulted Spots 

Occulted spots will appear in the transit light curve as 
a bump lasting roughly as long it takes the planet to move 
a distance 2(R p +R,) across the spot, assuming a circular 
spot with radius R,. For small spots (i?. << i? p ), this is 
roughly the transit ingress/egress time, or 6 minutes in 
the case of GJ1214b. The amplitude of a spot crossing 
event is determined by the fractional deviation in surface 
brightness occulted by the planet from the star's mean 
and could in principle be comparable to the transit depth 
itself for totally dark spots larger than the planet. 

We see evidence for spot occultations in the first two 
VLT transits shown in Fig. [4j in which the residuals from 



Cycle a 



T c (days) 1 



(days) 




9 
11 
21 
221 
228 
233 



2454966.525207 
2454980.748682 
2454983.909507 
2454999.713448 
2455315.794564 
2455326.857404 
2455334.759334 



0.000351 
0.000104 
0.000090 
0.000155 
0.000066 
0.000110 
0.000066 



a Time from the first MEarth 
transit, measured in units of the 
orbital period. 

b Mid-transit times T c are 
given as Barycentric Julian Dates 
in the Terrestrial Time system 
(BJDtdb)- To the best precision 
afforded by these transits, these 
times can be converted to a UTC 
based system via the approxi- 
mation BJDtjtc ~ BJDtdb — 
(32.184 + AO/86400 where N = 34 
is the number leap seconds be- 
tween UTC and TAT a ± the time 
of th ese transits (|Eastman et al.l 

[27lToI ). 

a smooth model show a 0.1% brightening at the start of 
the first transit and near the middle of the second tran- 
sit. These features last 5-10 minutes, as expected, and 
persist with comparable amplitudes in light curves gener- 
ated separately from the blue and red halves of the FORS 
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Fig. 7. — Deviations in the times of transit from the new best 
linear ephemeris, including transits from the V LT (filled black cir - 
cles), MEarth (open black circles), and work by Sada ct al. (2010) 
(open grey circles). 

spectra. Their presence is robust to choice of comparison 
stars. 

For the level of precision afforded by current data, 
we treat these possible spot occultations as excess red 
noise with correlations on the scale of the 6 minute 
spot crossing time. As such, the residual permutation 
method ( 34.2[) accounts for the uncertainty introduced 
by these features, regardless of their physical interpreta- 
tion. Spots that are not occulted will not be accounted 
for in t hese errors estimates, but we discuss their influ- 
ence in $6.2\ 

4.5. Unocculted Spots 

Unocculted spots will also have an effect on the plane- 
tary parameters. By diminishing the overall flux from the 
star while leaving the surface brightness along the transit 
chord unchanged, increasing the coverage of unocculted 
spots on the star will make transits deeper. This is op- 
posite the effect of occulted spots, which tend to fill in 
transits and make them shallower. Because we must ap- 
ply a nightly normalization to each light curve to avoid 
systematics and compare across telescopes, this depth 
change carries through to the implied 

To test whether unocculted spots are biasing individ- 
ual measurements of R p /R+, we repeat the fit and uncer- 
tainty estimation described in ^4. H and q4.2l but allow the 
8 transit epochs to have different values of R p /R±. As 
the visible spot coverage changes when the star rotates, 
we might detect changes in the inferred values of R p / R+. 
For this experiment we are interested in relative changes 
among the Rp/R* values, so we fix the other geometric 
parameters i?* and tu and the limb-darkening parame- 
ters to their best fit values to effectively collapse along 
those dominant degeneracies. 

The inferred i? p /-R* for each transit is shown in Fig. 
with uncertainties estimated from the residual permuta- 
tion method. Because the other parameters were fixed in 
this fit, the ensemble of points are free to move up and 
down slightly on this plot; the uncertainties shown are 
more relevant to comparisons between epochs. In ij6.2[ 
we will discuss our calculation of the predicted variations 
induced by unocculted spots consistent with the observed 
variability. 

5. LIMITS ON ADDITIONAL TRANSITING PLANETS 

After clipping out known transits of GJ1214b, we in- 
vestigate the light curve of GJ1214 for evidence of other 



Fig. 8. — Estimates of the apparent planet- to-star radius ratio 
at each epoch (black circles for MEarth/FLWO (open symbols) 
and VLT transits (filled symbols)). Via the residual permutation 
estimate, the error bars include the uncertainty due to possible 
presence of occulted spots. The predicted variation in the appar- 
ent Rp/Ri, due to the presence of unocculted spots (see i|6.2H is 
shown (grey), calculated directly from the nightly-binned MEarth 
photometry (points) or the best-fit sine curves to those data ( lines) . 

transiting planets in the system. The light curve contains 
3218 points spanning spring 2009 and 2010. 

To remove structured variability from the light curve 
before searching for transits, we employ an iterative fil- 
tering process_that^_combines a 2-day smoothed median 
filter (Aigrain fc Irwin! I2004D with a linear decorrelation 
against both external parameters (primarily the com- 
mon mode, meridian flip, seeing, and pixel position) 
and light curves of other field stars (|Kovacs et al.ll200"5t 
iTamuz et al.ll200l lOfir et al.l f2010h . Such filtering de- 
creases the light curve RMS from 9.6 to 4.7 millimagni- 
tudes per point. 

Using a variant of the box-fi tting least squares (BLS) 
algorithm ()Kovacs et al.l I2002D . we search the filtered 
light curve for periodic rectangular pulses over a grid 
of period and transit phase. At each test period, we fix 
the transit duration to that of a mid-latitude transit of 
a circularly orbiting planet for GJ1214's estimated stel- 
lar mass (0.157 M ) and radius (0.207 R Q ). With the 
typical MEarth precision and cadence, our sensitivity de- 
pends only weakly on transit duration; violations of these 
assumptions will not substantially penalize our detection 
efficiency. 

At every grid point, we determine the best fit transit 
depth D and the improvement A\ 2 over the D = null 
hypothesis, using a weighted least-squares method that 
includes red noise estimat ed from th e light curve itself 
using the oy formalism of iPont et al.l (|2006| ). Following 
iBurke et all (2006) we characterize candidates in terms 
of Ax 2 and / = max(Axf ach 

transit / A Xtotai)> which is the 
largest fraction that any one transit event contributes to 
the signal. 

The best candidate found in the clipped GJ1214 light 
curve exhibited a A% 2 of 31.3 and is shown in Fig. |9l 
To estimate the signific ance of this value, we employ 
the bootstrap method of I Jenkins et a l. (2002). Strictly 
speaking, the presence of remaining correlated noise in 
our filtered light curve means this method gives an over- 
estimate of the false alarm probability, but the compli- 
cated correlation structure of the light curve make a more 
accurate significance estimate difficult to calculate. In 
practice, we have found the Jenkins method to provide 
an appropriate limit, even for light curves with substan- 
tial red noise. We generate 10 3 (so we can estimate the 
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Fig. 9. — Left: The best candidate returned by the BLS search of the cleaned GJ1214 light curve. Bootstrap searches of simulated noise 
show the achieved Ax 2 =31.3 has a large 20% probability of false alarm. Right three: Examples of simulated transiting planets with periods 
near 10 days where the transits injected (black lines) into prefiltered light curves have crossed the detection threshold in the BLS search. 
All four panels show light curves phased to the period of the best fit BLS result (red dashed lines); the BLS does not necessarily recover 
the correct period. 



X 2 associated with a 1% false alarm probability from 10 
samples) transit-less light curves with Gaussian white 
noise and time sampling identical to that of the real light 
curve and performing the BLS search on these fake light 
curves. We find 20% of them show values of \ 2 > 31.3, 
suggesting our best candidate shown in Fig. [9] is not sig- 
nificant. 

These data place limits on the presence o f other tran- 
siting planets in the GJ1214 s ystem. Like iBurke et all 
ICroll et al.l (I2II0I and lBallard et all (|2010D .~^e 
simulate our sensitivity by injecting 8000 randomly 
phased, limb-darkened transits of 1.0, 1.5, 2.0, and 4.0 
R,0 planets with random periods. We then attempt to 
recover them with our transit search using objective de- 
tection criteria. To account for possible suppression from 
the filtering, we inject the transits into the raw light 
curves and reapply the filter at each iteration. We adopt 
two criteria for detection: A% 2 > 50 (corresponding to 
a formal probability of false alarm of 10 -5 by the boot- 
strap analysis) and / < 1.0 (to ensure at least two events 
contribute). 

To demonstrate this visually, we show in Fig. [5] ex- 
amples of injected transits for each relevant radius, ran- 
domly selected from among the simulated planets with 
injected periods of 10±1 days whose BLS results satisfied 
the detection criteria. Planets with periods near 10 days 
are of particular interest; given GJ1214's low luminosity, 
they would be in the star's habitable zone. Fig. [TU1 plots 
the fraction of injected planets that crossed the detection 
threshold as a function of period for each input planetary 
radius. The shape of the 4.0 R© curve in Fig. [Til] is driven 
largely by the / < 1.0 criterion requiring multiple events 
are observed, whereas the 1.5 R© curve is dominated by 
the need for sufficient in-transit S/N to get to A\ 2 > 50. 
With these criteria, our sensitivity falls below 50% be- 
yond periods of 15, 8, and 2 days for 4.0 R®, 2.0 R®, and 
1.5 R© transiting planets. We had no sensitivity to 1.0 
R© planets. 

Fig.QJJlalso shows the recovery fraction for 4.0 R© plan- 




Periods (days) 

Fig. 10. — The recovery fraction as a function of period for the 
labeled planetary radii (solid, colored lines) assuming the detection 
criteria A\ 2 > 50 and / < 1.0. For the 4.0 Rgj case, we also show 
the recovery fraction after lifting the / < 1.0 constraint (black 
dashed line); this is an estimate of the sensitivity to deep transits 
where only one event is necessary for a robust detection. 

ets without requirement that more than one event be ob- 
served. Neptune-sized transits of GJ1214 are so dramatic 
that they could be confidently identified by a single event 
(see Fig. [9]). Our simulations show that we are 90% sen- 
sitive to transiting Neptunes around GJ1214 out to 10 
days, and 80% sensitive out 20 days. For smaller planets 
requiring multiple events, however, a significant volume 
of parameter space remains unconstrained. 

6. DISCUSSION 

6.1. GJ1214 as a Spotted Star 

The slow rotation period we find implies a projected 
rotational velocity of 0.2 km s _1 that is well below the 
v sin i ~ 1 k m s" 1 detection limit of high-resolution rota- 
tion studies ([Browning et"aL 120101 : iDelfosse et al.l [19981 : 
iReinersI [20071 : iWest & Basrill2009l) but not inconsistent 
with lon g photometric period s detected for other field M 
dwarfs. iBenedict et al.l ([19981 ) estimated a rotation pe- 
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riod of 83 days for Pr oxima Cen, and recent photomet- 
ric work with A S AS d Kiraga fc Stepienl 120071) . H ATnet 



(|Hartman et all I2009T) . and MEarth (|Irwin et al.l 12010ft 

has confirmed the presence of many field M dwarfs with 
P rot > 10 days. 

Our P rot further implies that GJ1214 should ex- 
hibit signs of only weak magne tic activity (e.g. 
iReiners fc Basril l2007t iReinersI I2007D . Indeed, across 
three seasons of photometric monitoring, we see no evi- 
dence for flares in the MEarth bandpass, although their 
amplitude would be expected to be small in the near-IR. 
Activity induced chromospheric emission is not detected 
in either Ha or the Na I D doublet in the HARPS spec- 
tra used to measure the radial velocities presented in 
iCharbonneau et"aH (|2009f ). The relation between mag- 
netic activity and kinematic age ()West et al.l 120081) sug- 
gests that GJ1214 is > 3 Gyr old. 

We calculate GJ1214's ( U, V, W) space velocities 
(|Johnson fc So derblom 1987) in a left-handed system 
where U > in the direction of the Galactic anti-center 
to be (-47,-4,-40) km s -1 . These motions are consisten t 
with membership in the Galactic old disk (Leggctt 1992), 
lending further credence to an old age for GJ1214. 

6.2. Implications for Transmission Spectroscopy 

When inferring the transmission spectrum of a planet, 
one hopes to attribute changes in the transit depth across 
different wavelengths to atmospheric absorption by the 
limb of the planet. If the transmission spectrum is sen- 
sitive to 5 scale heights (H) of the planetary atmo- 
sphere, the amplitude of the transit depth variations 
are AD planct (A) = WHRp/B* or 0.001 if GJ1214 has 
an h ydrogen-rich atmosphere ( Millcr- Ricci fe Fortnevl 
2010). However, the presence of unocculted spots on 
the stellar surface can introduce transit depth variations 
AD 

spots (A) t) that are a function of both wavelength and 
time. To aid ongoing and future work to study GJ1214b's 
atmosphere, we use a simple model to estimate the am- 
plitude of the spot-induced contamination AZ? spots (A, t). 

We assume a fraction s(t) of the star's Earth-facing 
hemisphere is covered with spots; s(t) will change as the 
star rotates and the spots evolve. The observed out-of- 
transit spectrum -F 00 .t. (A, t) is a weighted average of the 
spectrum of the unspotted photosphere F Q {\) and that 
of the presumably cooler spotted surface F,{X): 



F . .t.(X,t) = [l-s(t)]F (X) + s(t)F.(X). 



(G) 



To simplify calculations, we neglect limb-darkening and 
treat each of the two components as having uniform sur- 
face brightnes s; tests with a limb-darkened spot model 
(|Dorrenl [1987) indicate only pathological cases could 
change the following results by more than 10%. 

When a planet with no atmosphere blocks light across 
a spot-free transit chord, it changes the relative weight 
of the two sources, causing the observed in-transit spec- 
trum Fi.t\\,i) to shift away from that of the unspotted 
photosphere toward that of the spots: 



^o.o.t.(A, t) 



= 1 - D(X,t) 



R, 



(7) 



I - ( ) - AD spots (X,t) 



1 - s(t) - 



F (X) + s(t)F.(X) 



[l-s(t)]F (X) + s(t)F.(X) 

Making the assumption that the total fraction of 
flux lost to the presence of spots is small (s(t)[l — 
F.(X)/F (X)] « 1), we solve to find 



AD spots {X,t) ps s(t) 



1 - 



F.(X) 



F (X) 



(8) 



If the quantity s(t)[l — F, (X)/F (A)] were not small, 
we probably would have observed larger amplitude and 
more frequent spot occultation events in the transit light 
curves in 2J We note the significant possibility that s(t) 
never reaches 0; that is, there may persist a population 
of symmetrically distributed spots that never rotates out 
of view. Such an unchanging population could cause us 
to overestimate the true value of R p /R+ by up to several 
percent. 

To match our observations of GJ1214's variability 
(AF„, .t,(A,t)/] 7 o.o.t.) to this model, we write s(t) — 
s + As(t) where s is the mean Earth-facing spot covering 
fraction and As(t) can be positive or negative. With the 
same assumption as above, we find 



AF n 



.(V) 



Fo.oA 



-As(t) 1 



F.(X) 

F (X) 



(9) 



In Sj3]we measured AF t (X,t) / F , o t . to have a peak- 
to-peak amplitude of 1% in the MEarth bandpass (715 < 
A < 1000 nm). Eq. [5] and [5] are equival ent to assuming a 
v alue of a = — 1 inlDesert et alj (|2011bD 's af\ formalism. 

iDesert et all (|201ibh show that measurements that 
rely on comparing photometric transit depths across mul- 
tiple transits could potentially mistake time-variability 
of A_D spo ts(i) for a feature in the transmission spectrum. 
From Eq. [5] and [91 we estimate the peak-to-peak time- 
variability of AD apota (t) to have an amplitude of 0.0001 
in MEarth wavelengths over the rotation period of the 
star. This spot-induced variability is comparable to, but 
smaller than, the estimated uncertainty from the tran- 
sits analyzed in this work and corresponds to an ap- 
parent change in planetary radius of AR p = 70 km or 
1/2 the sc ale height of an ^-dominate d atmosphere on 
GJ1214b (|Miller-Ricci fc Fortnevl [2Q10T l . In Fig. |we 
show the expected variation in the apparent planetary 
radius from unocculted spots and our individual 
measurements. 

Using blackbody spectra for F,(X) and F (X), we 
can extrapolate with Eq. [8] from the MEarth observa- 
tions (A w 0.85 /iin) to wavelengths accessible to Warm 
Spitzer. If we assume the spots are 300K cooler than the 
T c ff = 3000K stellar photosphere, we find AD(t) vari- 
ability amplitudes of 0.00004 and 0.00003 in Spitzer's 3.6 
and 4.5 /zm bandpasses. This allows robust comparison 
of transit depths between these wavelengths (Desert et 
al., in prep.). In the conservative limit that T, = OK, the 
variability amplitude would be achromatic. 

Even when comparison across different epochs is un- 
necessary, as is the case fo r spectroscopic observa- 
tions of individua l transits (e.g. lCharbonneau et al.|[2002l : 
iPont et all I2008D , unocculted spots can still introduce 
spurious wavelength features into the transmission spec- 
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trum. To place an upper limit on the amplitude of chro- 
matic changes in AZ? spot s(A) within a given transit, we 
imagine an extreme scenario where F m (X) is identical to 
-Fo(A) but with a very deep absorption line that does 
not appear in the unspotted spectrum, so the first fac- 
tor in Eq. |8j can vary with A between and s. Making 
the fairly conservative assumption that the population 
of spots that rotates in and out of view is comparable 
to the symmetric population (s < 4max[As]) we find 
that AL> spots (A) < 0.0003 for GJ1214 at wavelengths 
near 1 /xm. In practice, most features will show signifi- 
cantly lower amplitudes, although precise calculations of 
them and extrapolation to other wavelengths will require 
knowledge of the spot temperature and reliable model at- 
mospheres in the 2500-3000K temperature regime. 

Given the sensitivity of current generation instru- 
ments, the known population of spots on GJ1214 do 
not pose a significant problem for ongoing studies 
of GJ1214b's atmosphere. Future transmission spec- 
troscopy studies of GJ1214b, perhaps with the James 
Webb Space Telescope, comparing multiepoch, multi- 
wavelength transit dept hs and aiming to rea ch a precision 
of a D{\) = 0.0001 (see lDeming et afl 120091 ) will have to 
monitor and correct for the stellar variability. 

6.3. Limiting Uncertainties of GJ1214b 

As shown in Fig. [51 stellar spots currently play a very 
small role in limiting our understanding of the bulk mass 
and radius of GJ1214b. Here we address what other 
factors provide the limiting uncertainties in the planet's 
physical parameters. 

The 0.98 M ffi (15%) uncertainty on the planetary 
mass M p is the quadrature sum of 0.85 M ffi propagated 
from the measured radial velocity semiamplitude and 0.5 
Mq from the stellar mass uncertainty. Further spectro- 
scopic monitoring may reduce the former, but to fully 
reap the benefits of the radial velocities, the 12% er- 
ror on M+ should be improved. The current mass esti- 
mate is derived from an 2MASS photometry (2% uncer- 
tainty), an empirica l M^-mass relation (~ 10% scatter; 
iDelfosse et al.l I2000f) and the publishe d system parallax 
(77.2 ±5.4 mas: Ivan Altena et alll995[ ). Both improving 
Mif-mass relation for low-mass dwarfs and confirming 
GJ1214's parallax will be necessary to reduce GJ1214b's 
mass uncertainty. 

The 0.12 R.0 (5%) uncertainty on the planetary radius 
R p is already dominated by the 5% uncertainty in the 
stellar radius i?*, rather than the light curve parame- 
ters. This is currently constrained by the combination 
of the stellar densit y p+ that is measured directly from 
transit light curves ()Seager fc Mallen-Orn elas 2003|) and 

the estimated mass. Even though cx (M^/p*) 1 / 3 , and 
so is relatively insensitive to uncertainty in M*, if we fix 
M* = 0.157 ± 0.019M Q we find the errors on R± and 
R p shrink by a factor of two. Improving the estimate 
of GJ1214's mass is also the best way to improve our 
measurement of the radius of the planet. 

We reiterate here that our measured out-of-transit pho- 
tometric modulation probes only the spatially asymmet- 
ric component of the stellar spot distribution rotating 
around the star. If the star hosts a subtantial, un- 
changing, spatially symmetric population of unocculted 
spots, it will bias estimates of the true planetary radius 



(jCzesla et al. 2009) too high. If the symmetric and asym- 
metric components of the spots are comparable, such a 
bias will be at the percent level of in R p , smaller than 
the current uncertainties. 

6.4. Metallicity of GJ1214 

Several authors have recently developed empiri- 
cal photometric calibrations to estimate M dwarf 
metalliciti es from abso l ute Mk magnitude and V 



Johnson fc Appsl [2009; 



K color dBonfils et all |200~ 

iSchlaufman fc Laughlinl [2010D . As exoplanet surveys 
like MEarth lavish more attention on M dwarfs as ex- 
oplanet hosts, such studies hope to address whether 
the giant planet vs. stellar metallicity correlation 
seen for AFGK ste llar hosts (|Fischer & Valentil 120051 : 
Uohnson et al.l 
smaller stars. 

tion does not persist down to Neptune-mass planets, but 
more data are needed. 

Our improved estimate of V — 14.71 ± 0.03 differs 
si gnificantly from the V = 1 5.1 central value published 



20101) extends to smaller planets and 



Sou sa et al.l ()2008l ) suggest the correla- 



Charbonneau et ail (|2009t ) but was the value used in 



the most recent a nalysis of this photome t ric me tallic- 
ity calibration by ISchlaufman fc Laughlinl ([2010) , who 
found [Fe/H] = +0.28 fo r GJ1214. This an a lysis a grees 
quite well with work by iRoias-Avala et all ([2010) that 
estimates empirically calibrated metallicities from alkali 
metal lines in moderate resolution if-band spectra and 
finds [Fe/H] = +0.39 ±0.15 for GJ1214. This lends in- 
cremental evidence towards the persistence of the mass- 
metallicity correlation down to super-Earths around M 
dwarf hosts. 

7. CONCLUSIONS 

We have measured long-term photometric variability 
on GJ1214 to have a 1% peak-to-peak amplitude in the 
MEarth bandpass (715-1000 nm) and a long rotation pe- 
riod, most likely an integer multiple of 53 days. Fitting 
very high precision light curves from the VLT, we find 
likely instances of GJ1214b crossing small spots during 
transit. Treating these occultation events as correlated 
noise, we find parameters for the planetary system that 
are consistent with previous work. 

We estimate the amplitude of time- variable changes in 
the apparent radius of the planet due to the observed stel- 
lar variability as AD spo t s (t) = 0.0001 and place an upper 
limit of AD spots (A) < 0.0003 on possible spot-induced 
spectral features in the planet's transmissio n spectrum- 
Stellar spots do not limit current studies (e.g. lBean et al.l 
120101 : iDisert et aLll2011a|) . but could be important for fu- 
ture studies of GJ1214b with JWST. 

Using two years of MEarth data, we have placed lim- 
its on the presence of other transiting planets around 
GJ1214. With 90% confidence, we rule out the presence 
of Neptune-radius transiting planets in orbits shorter 
than 10 days but cannot place strong constraints on 
planets smaller than 2.0 R® at such long periods. In 
a system where a 1.0M ffi planet in a 2:1 mean mo- 
tion resonance would create 100 second perturbations to 
GJ1214b's transit times, we find no evidence for tran- 
sit timing variations larger than 15 seconds. Further 
searches of the GJ1214 system for potentially habitable 
planets smaller and cooler than GJ1214b continue to be 
warranted. 
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